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Abstract 

We study the blow up solutions of a semilinear reaction diffusion system cou- 
pled in both equations and boundary conditions. The main purpose is to un- 
derstand how the reaction terms and the absorption terms affect the blow-up 
properties. We derive the lower and upper bound for the blow-up rate, and find 
the blow-up set under certain assumptions. 



1 Introduction 

In this paper, we consider the following parabolic system 

u t = Au + X 1 e v , v t = Av + X 2 e u , (x,t) G B R x (0,T), 

% = e u , (x,t)edB R x(0,T), 
u(x,0) = u (x), v(x,0) = v (x), x G Br, 

where Ai,A2 > 0, uo,vo are nonnegative, radial nondecr easing, smooth functions 
and satisfy the conditions 

if = e *°> W = e " ' x€dB R , ] 

Auo + e Vo > 0, Av + e u ° > 0, x G B R , > (1.2) 

u 0r (\x\) > 0, v 0r (\x\) > 0, x G Br. ) 

The problems of semilinear systems coupled in both equations and boundary 
conditions have been studied very extensively over past years in case the reaction 
terms and boundary conditions are of power type functions, for instance in [3], it 
was considered the solutions of the following system 



u t = u xx + vP\ vt = v xx + U P*, (x,t) G (0, 1) x (0, T), 

u x (l,t) = vi\ v x (l,t) = u^, iG(0,T), 

u x (0,t)=0, v x (0,t)=0, iG(0,T), 

u(x, 0) = uq(x), v (x, 0) = vq(x), i£ [0,1], 



;i.3) 
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where pi,P2, 9i> 92 > 0, and uo, vq are radial nondecr easing, positive smooth func- 
tions satisfying the conditions 

«0x(0) = v Ox (0) = 0, u 0x (l) = <(1), v 0x (l) = u q 2 (l). 

It was shown that if 

max{p 1 p 2 ,p 1 q 2 ,P2qi,qiq2} < h 

then the solutions of problem (|1 .3H exists globally, otherwise every solution blows 
up in finite time. Moreover, the blow-up occurs only at x = 1 and the blow-up 
rate estimates take the following form 

Ci(T-t)- Q <u(l,t) <C 2 (T-t)- a , t € (0,T), 

C 3 (T-t)-P <v(l,t) <C 4 (T-t)-P, te(0,T), 

where 

a = a(pi,p 2 ,qi,q 2 ), /3 = £(pi,P2, Qi, 

In [9], it was considered the critical exponents for a system of heat equa- 
tions with inner absorption reaction terms and coupled boundary conditions of 
exponential type, namely 

ut = An - a ie P lu , v t = Av - a 2 e P2V , (z, t) G fi X (0, T), } 
%=^\ % = z q2U , (x,t)edSlx(p,T), \ (1.4) 

u(x,0) = uq(x), v(x,0) = Vq(x), X E U, ) 

where £1 is a bounded domain with smooth boundary, pi,p 2 > 0, gj,Oj > 0,i = 
1,2, uo,t>o are nonnegative functions that satisfy 

or/ 

It was shown that if 

1/n > 0, or l/r2 > 0, 

where 

9192 - iPlP2 ' 9192 - 2P1P2 ' 

then the solutions of problem (jl.4p with large initial data blow up in finite time. 

The main purpose of this paper is to derive the upper and lower blow-up rate 
estimates for problem (jl.ip and to study the blow-up set under some restricted 
assumptions. 
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2 Preliminaries 



Since the system is uniformly parabolic, also the reaction and the boundary 
conditions terms are smooth functions and the initial data satisfy the compat- 
ibility conditions, therefore, the local existence and uniqueness of the classical 
solutions of problem are known by standard parabolic theory (see [3]). On 
the other hand, for any initial data (uo,vo), the solution of this system has to blow 
up in finite time and the blow-up set contains the boundary (8Br), and that due 
to the comparison principle [7j and the known blow-up results of problem (jl.ip . 
where Ai = A2 = 0, which has been studied in [6]. 

The next lemma shows the properties of the classical solutions of problem 
(11.11) . We denote for simplicity u(r, t) = u(x,t). 

Lemma 2.1. Let (u,v) be a classical solution to problem U.l\) . Then 

(i) (u,v) is radial and u, v > in Br x (0,T). 

(ii) u r ,v r > in [0,22] x (0,T). 

(iii) ut,iH> 0, in B R x (0,T). 

Next, we prove the following lemma, which shows the relation between u and 

v. 

Lemma 2.2. Let (u,v) be a solution to problem M.l}) , there exist M > 1 such 
that 

e v <Me u , e u <Me v , (x,t) € Br x [0,T). (2.1) 

Proof. Let 

J(x,t) = Me"^ - e v ^\ (x,t) € Br x (0, T), r = |x|. 
A direct calculation shows 

J t = Me u u t -e v v t , 

J r = Me u u r -e v v r , (2.2) 
J rr = Me u u rr + Me u u 2 r - e v v rr - e v v 2 



r ■ 



Thus 



J t - J rr - - — - J r = Me u u t - e v v t - Me u u„ - Me u ul + e v v rr + eV 
r 

—Me u u r H eV 

?2 — 1 7Z — 1 

= Me u ktj — u rr u r ] — e v [vt — v rr v r ] 

r r 

-Me u u 2 r + e v v 2 r 
= Me u [X l e v ] - e v [X 2 e u ] - Me u u 2 r + e v v 2 . 
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Prom it follows that 



Ur = A^ [VreV + Jr] > 

U r = Jp^I^r e2V + 2eVv r J r + J r\- 

Therefore, 

J t -AJ= (AiM - A 2 )e ,t+t ' + [e 1 ' - -^R 2 - \-^-v r + -r^ J r ] Jr- 
v 7 L Me ul r v Me u Me u 1 

Clearly, 

e 2v T 



Me" Me" 
Therefore, the last equation can be rewritten as follows: 

J t -AJ- bJ r -cJ= (AiM - A 2 )e" +1 ' > 0, (a?, t)e B R x (0, T) 

provided M > A2/A1, where, 

^e" 1 T , e v 2 

= — I TT v r + TT Jrh c = TT v r ■ 

l Me u Me u 1 Me u r 
It clear that, b, c are continuous functions and c is bounded in Br x (0, T*), for 

Moreover, 

<9J 

^■Ue9Bfl = [Me u u r - e v v r ] 

= Me u+V - e u+v = [M- l]e u+v > 0, 

and 

J(x, 0) = Me"° - e V0 > 0, x £ ~Br 

provided M is large enough. 

From above and Proposition [7j, it follows that 

J>0, in B R x[0,T). 

Similarly, we can show that the function H = Me" — e u is nonnegative in Br x 
[0,T). □ 
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3 Blow-up Rate Estimates 

In this section we consider the upper and lower blow-up rate estimates of solutions 
for problem (fTTTj) with (fOj) . 

Theorem 3.1. Let u be a blow-up solution solution of problem j fi.ip 
X\ = A 2 = A, T is the blow-up time. Assume that uq,vq satisfy 

«or(r) - ^e" o(r) > 0, %(r)-^e u »W>0, re[0,E]. (3.1) 

T/ien there is a positive constant c such that 

logc--log(T-t)<u(R,t), \ogc--log(T-t)<v(R,t), te{0,T). 

Proof. Define the functions J±, J 2 as follows: 

Jl(x,t) = U r (r,t) - L e <r,t)^ j^ x t) = _ J_ e u(r,t)_ 

K K 

A direct calculation shows 

r n — 1 

Jit = u r t ~ -^e v [v rr -\ v r + \e u ], 

R r 



T 

v 

n — 1 n — 1 

J\ rr = [u r t U rr H 5 — U r — Xe V V T ] 



r 2 
--[e v v rr + e v v 2 r ] - jf v v r . 



From above it follows that 

n — 1 . n — 1 



- Jirr — J ir = — K - -^e"] + Ae"[u r - — e u ] + ~^e v v; + — e*V 



Thus 

Jit - AJi + ^ Ji - Ae 1 ' J 2 = ^e^ r 2 + \e v v r > 0, 

R R 

for (x, t) e B R x (0, T) n {r > 0}. 
In the same way we can show that 

n — 1 

J 2i - A J 2 + —g- J 2 - Ae u J x > 0, (x, t)eB R x (0, T) n {r > 0}. 

Clearly, from (|3.ip . it follows that 

Ji(x,0), J 2 (x,0) > x € Br. 
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And 

Ji(0,t) = Ur(0,t) > 0, J 2 (0,t) = v r {0,t) > 0, 
Ji(R,t) = J 2 (R,t) = 0, te(0,T). 

Since, the supremums of the functions Ae M , Ae^ and -f^ (on x (0, t] for t < T) 
are finite, therefore, from above and maximum principle, it follows 



Moreover, 



This means 



Thus 



Ji,J 2 >0, (x,t) eB R x (0,T). 



r 1 

(«rr - ~ -j^ eV )\dB R < 0. 

71 — 1 T 

ut < ( u r + Xe v + -e v v r + -^la^. 

r ti ti 



Which implies that 

Ut (R it ) < V^_± e v(R,t) + Ae ,(R,t) + e v(R,t)+u(R,t) + 1 t G (0jT) 

From the last inequality and Lemma 12.21 it follows 

Ut(jR) t) < ^Me^-*' + AMe^" + Me 2 ^ 4 ' + ^ e «(*.*), f g (0, T). 

Thus, there exist a constant C such that 

u t (R,t) <Ce 2 < R *\ t€(0,T). 

Integrate this inequality from t to T and since u blows up at R, it follows 

£_^< e «(fl,t) te(0,T) 

(T-t)3 " 

or 

logc--log(T-i)<u(i2,<), tG(0,T). 
We can show in a similar way that 

logc-hog{T-t)<v(R,t), te(0,T). 



□ 



Next, we consider the upper bounds 

Theorem 3.2. Let u be a blow-up solution solution of problem T is 

the blow-up time. Then there is a positive constant C such that 

u(R,t) <logC-log(T-i), v(R, t) < log C — log (T — t), te(0,T). 

Proof. Define 

M(t) = m&xu(x,t), N(t) = m&xv(x,t). 

Br Br 

M(t),N(t) are increasing in (0,T) due to the 

ut,v t >0, (x,t) eB R x (0,T). 

For < z < t < T,x £ Br, as in [1], the integral equation for problem (jl.ip with 
respect to u can be written as follows 

u{x,t) = [ T(x-y,t- z)u{y,z)dy + \ 1 [ [ T(x - y,t - T)e v(y ' r) dydr 

J Br Jz JBr 

+ 11 T{x-y,t-T)e v ^ds y dT 

Jz JSr 

-[ [ u (y, T )-7r~(x - V,t - T )ds y dT, 
Jz Js R drjy 

where T is the fundamental solution of the heat equation, which takes the form: 

1 Id 2 
r(x, t) = -r^r exp[— !-M. (3.2) 

Letting x — > 8Br and using the jump relation, [2], for the fourth term on the 
right hand side of the last equation, we obtain 



1 / \ 
-u(x,t) 



f F(x-y,t-z)u(y,z)dy + X 1 f f T(x - y,t - r)e^ r) dydr 
+ [ [ r(x-y,t-T)e v{v ' T) ds y dT 

Jz JSr 

f* f dT 

~ / u(y,T)-—(x-y,t-T)dsydT, 
Jz Js R dr] y 



'Sr 

for x G dB R , < z < t < T. 

Since u, v are positive and radial, it follows 



T(x- y,t- z)u(y,z)dy > 0, 

Br 

t r rt 



f f e v(v,r) T ( x -y t t- T)ds y d T = f e v ^ [f r(x-y,t- T)ds y ]dT. 

J z J Sr J z J Sr 



Thus 



\M{t) > I e N ^[f r(x-y,t-T)ds y ]dT 

2 Jz JSr 

r l r dr 

- M(t)[ \ — {x-y,t-r)\ds y }dT, x £ Sr, < z < t < T. 
Jz Js R ar ly 

It is known that (see [2]) for < t 2 < t 2 , these is C* > such that 

dT C* 1 

W x ~ y,h ~ h)l ~ W^W ' \x- y \(n + ^-«y x > y G s ">° € ^ ^ 

Choose 1 — ^ < fi < 1, from [2], there exist Ci > such that 

J Sr \x-y\("+±-^-<r) V 
Also, if close to t 2 , then there exist a constant c such that 

/ T(x -y,t 2 - h)ds y > 



Thus 

1 r* e N ^ r* m(t) 

-M(t) > c / ^=dr - C / , [ dr. 

2 y, ~ Jz Vt^T h \t - r|M 

Since for < z < r < t < T, it is clear that M(r) < M(i), thus 
1 /"* e^W 

-M(t) > c J ^—dr - C{M{t)\T - z\ l ~». (3.3) 
Taking z so that C{\T - = 1/2, it follows 

M(t) > c jf -^^dr = A(t). (3.4) 

Clearly, 

A(t) = c W=t 

Prom Lemma [2.2| there exist a constant k > 1 such that the last equation becomes 
which leads to 

> 



e A ~ J t k y/T- 



Clearly, 



ft e N(r) rt e N(r) 

A(T) = lim c / dr = c lim dr = oo 

t^r J z J z t—>t Vr^T 



This leads to 

1 2c 



> — Vr^t. 



Therefore, there exist a constant Co > such that 



e A(t) < Co = z<f<T (35) 



On the other hand, for to = 2t — T (Assuming that t is close to T), 

-t p N{r) 
't 

Combining the last inequality with (|3.5p . yields 

C 



rt N(t) rt 1 

A(t) > c / -. dr > C e N ^ / dr = e N{t °hc{V2 - \)y/T=t. 



which leads to 



> e N{ta) 2c{V2- l)VT=t, 
Co 



e N(t ) < 

" c(V2 - 1)(T - t )' 

Thus there exist a constant C such that 

e *(*)< <t<T 
" (T-t)' 

or 

v(R, t) < log C — log (T — t), te(0,T). 
In the same way we can show 

u(R,t) < log C- log (T-i), ie(0,T). 



□ 



Remark 3.3. From Theorems 13. Hl3T2| we conclude that, the upper blow-up rate 
estimates of problem (jl.ip are coincident with the upper blow-up rate estimates 
of the zero Dirichlet problem for the semilinear system in (jl.ip . while the lower 
blow-up rate estimates of problems (|l.ip are coincident with the lower blow-up 
rate estimates of problem (jl.ip . where Ai = A2 = (see [1]). 
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4 Blow-up Set 



We consider next the blow-up set for problem (jl.ip . under some restricted as- 
sumptions on Ai, A2. 

Theorem 4.1. Let (u,v) be a blow-up solution to problem M.l}) . Assume that 
the following condition is satisfied 



A[4tf (n+D + 1] < min j-, [R2 + ^ -2— - + ^ ^ e 

(4.1) 

where T is the blow-up time, C is given in Theorem \3.SX A = max{Ai, A2}. Then 
there exist a positive constant A such that 

u(x, t) < log[ 1 — ^2], v(x, t) < log ' 



-A(R 2 -r 2 ) 21 ' v ' ; " bl A{R 2 -r 2 ) 2J ' 
/or efifix (0,T), 
Proof. Define the functions Z\ , z 2 as follows 

zi(x,t) = z 2 (x,t) = log [Av{x) + B{T _ t)] , (x,t) e Sfl x (0,T), (4.2) 



where u(x) = (i? 2 — r 2 ) 2 , r=|x|, > 0, ^4 > A. 
A direct calculation shows: 

z it - Azi - Xie Z2 > z lt - Azx - Ae Z2 > 0, in B R x (0, T), 
z 2t ~ Az 2 - X 2 e Zl > z 2t - Az 2 - Ae Zl > 0, in B R x (0, T) 

provided 

B > A[4R 2 {n + 1) + 1]. 

Moreover, 

zi(x,0) = log [Av(x ] +BT] > log [AR i\ BT] > u(x,0), xeB R , 
z 2 (x,0) = log [Av{x ] +Bn > log [ARi \ BT] > v(x,0), xeB R 



C" i? 2 + 4(n + 1)T , R 2 + 4(n + l)T 

From (1131) and Theorem E21 it follows that 

zi(i?, t) > t), z 2 (R, t) > v(R, t), t€ (0, T), 
zi(a;, 0) > u(x, 0), Z2(x, 0) > v(x, 0), x £ B R . 



(4.3) 



(4.4) 



and 

zi(i2,t)=2 2 (12,t)=log 5 ^>logpg i j, ie(0,T) (4.5) 

provided 

, 1 4(n + l) I, 4(n + l) 

B < mini — , 



(4.6) 
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Prom (|4.3p . (|4.6p and the comparison principle jS], it follows that 
> z 2 0M) > v(a;,i), (i,t)eB s x(0,T). 

Moreover, from ()4.2f) 




(4.7) 



□ 



Remark 4.2. From (|4.7p . we conclude that, for problem (jl.ip with (|4. 1 1) . any 
point x € Br cannot be a blow-up point, therefore, the blow-up occurs only at 
the boundary. This means, if Ai, A2 are small enough, then the blow-up set is the 
same as that of (jl.ip . where Ai = A2 = (see [6]). 
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